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ABSTRACT 

Gravitational radiation tends to drive gravity modes in rotating neutron stars unstable. 
For an inviscid star, the instability sets in when the rotation frequency is about 0.7 
times the corresponding mode frequency of the nonrotating star. Neutron stars with 
spin frequencies ^ 100 Hz are susceptible to this instability, with growth time of order 
years. However, it is likely that viscous dissipation suppresses the instability except 
for a narrow range of temperatures around 10^ K. We also show that the viscosity 
driven instability of g-modes is absent. 
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1 INTRODUCTION 

It is well known that globally nonaxisymmetric instabili- 
ties can develop in rapidly rotating neutron stars. Of par- 
ticular interest is the secular instability driven by gravita- 
tional radiation, as first discovered by Chandrasekhar (1970) 
for the bar-mode (m = 2) of an incompressible Maclaurin 
spheroid, and later shown by Friedman & Schutz (1978a,b) 
to be a generic feature of rotating stars. The Chandrasekhar- 
Friedman-Schutz instability (hereafter CFS instability) oc- 
curs whenever a backward-going mode (with respect to the 
rotation) in the corotating frame of the star is "dragged" 
by rotation to become forward-going in the inertial frame. 
A CFS unstable mode has negative energy (in the inertial 
frame), and gravitational radiation makes the energy even 
more negative. 

The CFS instability has been extensively studied for 
f-modes since it is widely believed that the maximum angu- 
lar velocities of neutron stars are determined by instabilities 
of these modes (e.g., Ipser & Lindblom 1990, 1991; Cut- 
ler & Lindblom 1992; Lindblom 1995; Yoshida & Eriguchi 
1995; Stergioulas & Friedman 1998). The f-mode instabili- 
ties occur at rotation frequencies comparable to the max- 
imum "break-up" frequency (~ 1000 Hertz). It was found 
that viscous dissipation in the neutron star stablizes almost 
all f-modes except in certain limited temperature regimes. 
Recently Andersson (1998) pointed out that r-modes (pre- 
viously studied by Papaloizou & Pringle 1978, Provost et 
al. 1981, Saio 1982, etc.) are also destablized by gravita- 
tional radiation (see also Friedman & Morsink 1998). More- 
over, the couplings between gravitational radiation and r- 
modes (mainly through current multipole moments) are so 
strong that viscous forces present in hot young neutron stars 
are not sufficient to suppress the instability (Lindblom et 



al. 1998; Andersson et al. 1998a). It is therefore likely that 
gravitational radiation can carry away most of the angular 
momentum of a nascent neutron star through the unsta- 
ble r-modes, and the gravitational waves generated in this 
process could potentially be detectable by LIGO (Owen et 
al. 1998). The implications of the r-mode instability for ac- 
creting neutron stars have also been discussed (Andersson 
et al. 1998b; Bildsten 1998; Levin 1998). 

A neutron star also possesses p-modes and g-modes, 
whose stability properties have not been studied. The p- 
modes are similar to the f-modes, except that, as compared 
to the f-modes, they have higher frequencies, weaker cou- 
plings to the gravitational radiation and stronger viscous 
dampings. Thus the p-mode instabilities, if exist, are cer- 
tainly less important than the f-mode instabilities. The sit- 
uation with g-modes is not clear, however. Gravity modes 
in neutron stars arise from composition (proton to neutron 
ratio) gradient in the stellar core (Reisenegger & Goldreich 
1992; Lai 1994), density discontinuities in the crust (Finn 
1987; Strohmayer 1993) as well as thermal buoyancy asso- 
ciated with finite temperatures, either due to internal heat 
(McDermott, Van Horn & Hansen 1988) or due to accre- 
tion (McDermott & Taam 1987; Bildsten & Cutler 1995; 
Strohmayer & Lee 1996; Bildsten & Gumming 1998). For 
a nonrotating neutron star, the typical frequencies of low- 
order g-modes are around 100 Hz and are smaller for higher- 
order modes. Such low frequencies imply that, contrary to 
f-modes, g-modes may be destablized by gravitational wave 
when the star rotates at a rate much smaller than the break- 
up spin frequency. Moreover, since g-modes possess small, 
but nonzero mass quadrupole moments, it is not clear a 
priori whether the driving due to gravitational radiation is 
stronger or weaker for g-modes than that for r-modes (which 
have zero mass quadrupole moment). The purpose of this 
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paper is to address these issues related to the secular insta- 
bility of g-modes. 

Our paper is organized as follows. We start in Section 2 
by summarizing the essential properties of g-modes in nonro- 
tating neutron stars and defining notations. In Section 3 we 
calculate the effect of rotation on the mode frequencies and 
determine the onset of CFS instabilities of g-modes for invis- 
cid stars. In Section 4 we evaluate the driving/damping rates 
of the g-modes due to gravitational radiation and viscosi- 
ties. Section 5 considers the possibility of viscosity-driven 
g-mode instability. In Section 6 we compare the instabilities 
of g-modes and other modes (f and r-modes) and discuss 
some possible astrophysical applications of our results. 



2 PREPARATIONS: CORE G-MODES OF 
NEUTRON STARS 

In this paper, we shall focus on the core g-modes, with 
buoyancy provided by the gradient of proton to neutron ra- 
tio in the neutron star interior (Reisenegger & Goldreich 
1992). These modes exist even in zero-temperature stars. 
The other types of g-modes (due to crustal density discon- 
tinuities or thermal entropy gradient) either have weaker 
restoring forces or are confined to the outer region of the 
star and therefore are less susceptible to the CFS instabili- 
ties as compared to the core g-modes (see below) . 

We construct equilibrium neutron star models by solv- 
ing Newtonian hydrostatic equations. Centrifugal distortion 
to the stellar structure is neglected as we shall only con- 
sider stars with rotation rates far below the break-up value. 
The equations of state (EOS) used for our neutron star 
models and g-mode calculations have been described in Lai 
(1994). The nuclear symmetry energy plays an important 
role in determining the proton fraction x = Up/n (where n 
is the baryon number density and Tip in the proton num- 
ber density) and hence the frequencies of g-modes. Stable 
stratification of the liquid core (thus the existence of core g- 
modes) arises from the a;-gradient, and from the fact the f3- 
equilibration time is much longer than the dynamical times 
of the modes. The Brunt- Vaisala (angular) frequency A'' is 
given by 

9 fdp^ 



N 



'dx\ 

p \dxjp \dr J 



(1) 



where <; > is the local gravitational accelaration. As 
an estimate of the mode frequency, we neglect the nu- 



nuc ~ 2.8 X 10^* g cm ^ is the nuclear density. The Brunt- 



cleon interaction, and find x 

p. 

Vaisala frequency is then 

1/2 



5.6 X 10 (p/pnuc), where 



Vlop" 



(2) 



Thus the core g-modes typically have frequencies less than 
100 - 200 Hz. 

In the following sections, we shall consider two repre- 
sentative models based on a parametrized EOS derived from 
Wiringa, Fiks & Fabrocini (1988), and labeled as UU and 
AU in Lai (1994). Figure 1 depicts the profiles of the Brunt- 
Vaisala frequency, I'm ~ N/{2tt), as a function of radius 
inside the star for the two models. We see that z^jv = 
at the center (where g = 0) and increases outward be- 
cause of decreasing sound speed Cs and increasing g (de- 
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Figure 1. Profiles of the Brunt- Vaisala frequency vjv (heavy 
curves) and the Lamb frequency vj^ (light curves) for two neu- 
tron star models: UU (solid curves) and AU (dashed curves). We 
have set vjsi = outside the core (p ^ 10^"* g cm~^). The dips 
in the vj^ curves at rjR ~ 0.94 correspond to neutron drip at 
P - 



: 4 X lO" g cm-3. 



spite that x decreases outward). We have set vm = out- 
side the core (p ^ 10^* g cm~"^)^ and thus eliminated 
the crustal modes (see below). For completeness. Figure 1 
also shows the Lamb frequency (or the acoustic cut-off fre- 
quency), PL ~ + 1) Cs/ (2-Kr), with I = 2. Gravity wave 
(with frequency i^) can propagate only in the region satisfy- 
ing u < vl and v < vn- For a typical v ~ 100 Hz, the wave 
propagation region lies in the range 0.5 ^ r/R 0.9. For 
convenience, the basic properties of the modes for nonrotat- 
ing stars are summarized in Table 1. 

The crustal discontinuity modes have angular frequen- 
cies approximately given by u ~ {g / R)^^"^ {/S.p / p)^^'^ . With 
typical density discontinuity ^p/ p ~ 1%, this gives a fre- 
quency of order 200 Hz (Finn 1987). Although these modes 
can have frequencies greater than the core g-modes, the 
crustal modes are confined in the outer region of the star and 
therefore have very small mass quadrupole moments^. The 
thermally driven g-modes (McDermott et al. 1988) have fre- 
quencies less than 50 Hz for T <, 10^" K, indicating smaller 
restoring forces as compared to the core g-modes. We will 
not consider the crustal modes and thermal modes in the 
rest of this paper. 



* The transition from the core (where only free nucleons and 
electrons exist) to the inner crust (where heavy nuclei co-exist 
with free neutrons) is somewhat ambiguous. The inner crust is 
described by the equation of state of Baym, Bethe & Pethick 
(1971) while the core is described by that of Wiringa et al. (1988). 
We have found the g-mode properties to be insensitive to the 
actual transition radius within the allowed range, 
t This is indicated by the fact that the mode damping time due 
to gravitational radiation is rather long, exceeding 10'* — 10^ years. 
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Table 1. Core g-modes of nonrotating neutron stars 



Model 


M (Mq) 


R (km) 


Pc (g cm 3) 


mode ^ 


u (Hz) 


15^22 1 [ 


1 


UU 


1.4 


13.48 


6.95 X 10" 


gi 


148 


6.2 X 10" 


-4 










g2 


97 


1.8 X IQ- 


-5 


AU 


1.4 


12.37 


8.18 X 10" 


gl 


72 


1.1 X IQ- 


-4 










g2 


52 


9.9 X 10^ 


-6 



" The subscript n in gn specifies the radial order of the mode. The angular order is Z = 2. 
^ Note that the values of |<5D22| quoted here differ from those in Lai (1994) by up to 80%. This is because we 
adopt Cowling approximation in this paper while in Lai (1994) the full equations are solved; see discussion 
at the end of §3. 



3 ONSET OF INSTABILITY: INVISCID STARS 

For neutron stars consisting of ideal fluid (neglecting viscos- 
ity), it is possible to determine the onset of CFS instability of 
g-modes without calculating the growth/damping rates due 
to dissipative processes. Consider a mode with azimuthal an- 
gular dependence e*™*^ and time dependence e"^''^, where uj,. 
(> 0) is the mode angular frequency in the rotating frame. 
The rotation angular frequency is f2s (> 0). A mode with 
m > propagates in the direction opposite to the rotation, 
with pattern speed {—LOr/m) (in the rotating frame). In the 
inertial frame, the mode frequency is 



(3) 



When uJi < 0, the mode attains a prograde pattern speed, 
— {oJi/m), in the inertial frame, and the CFS instability sets 
in (Friedman & Schutz 1978a,b). 

Since we shall consider Qs comparable to Ur, a per- 
turbative treatment of the rotational effect is inadequate 
(see the last paragraph of this section). We shall adopt the 
so-called "traditional approximation" (Chapman & Lindzen 
1970; Unno et al. 1989; Bildsten et al. 1996), where the Cori- 
olis forces associated with the horizontal component of the 
spin angular velocity are neglected. This allows the mode 
eigenfunction to be separated into angular and radial com- 
ponents. Note that the traditional approximation is strictly 
valid only for high-order g-modes for which ujr N and 
3> (where ~ ~ the horizontal Lagrangian 

displacement and is the radial displacement), and for 
fls <^ N (so that the radial component of the Coriolis force 
can be neglected compared to the buoyancy force). While 
these conditions are not valid everywhere inside the star 
(e.g., near the center is of the same order as ^h, and 
uir,Qs > N; also, near the surface, uir,Qs > N), in the re- 
gion (0.6 ^ r/R 0.9) where the mode propagates, they 
are reasonably satisfied (see below). 

Neglecting the perturbation in gravitational potential 
(Cowling approximation), the radial Lagrangian displace- 
ment and Eulerian pressure perturbation can be written as 



5P{t) = 5P{r)Hj^{9)e''^t, (4) 



where Hjm{0) is the Hough function (j'm are the indices 
analogous to Im in the spherical harmonic Yjm), satisfying 
the Laplace tidal equation (Chapman & Lindzen 1970; Bild- 
sten et al. 1996): 



d 



m Hi 



(l-M2)(l-gV) 



gm(l + q^fjp)Hj 



(5) 



(1-gV^)^ 

with /i — cos 6 and q = 2Qs/i^r- The eigenvalue. A, depends 
on m and q. For q ^ (the nonrotating case), the function 
Hjm{6)e^""^ becomes Yim{6, <j)) while A degenerates into + 
1). The other eigenfunctions of the mode are related to 5P 
and via 



5p{v) = 5p{r)H,^e-'' = ( ? + —6 

Cs g 



?0(r) = 



dH 



de 

dH. 



H : — -lii 



sin 6 



( oHjrn , m 

I'^'^^ + ihT?^- 



,(6) 
(7) 
(8) 



where Cs is the sound speed and £,±{r) = 5P{r) / {pruj^). 
Separating out the angular dependence, the fluid continuity 
equation and Euler equation (in the rotating frame) reduce 
to a set of coupled radial equations: 

d 



dr 



-SP 



rSP + piiut - N')^r 



^ir'^r) = ^{r%)-^6P+^5P 
dr d dp u^p 



(9) 



(10) 



The properties of the Hough function have been exten- 
sively studied (Longuet-Higgins 1967). We have adopted the 
numerical approach of Bildsten et al. (1996) in our calcula- 
tion. For given q and m, the eigenvalue A is obtained by 
solving equation (jsj). To obtain the eigen- frequency Ur, we 
solve the radial equations (P)-(p^ subjected to the follow- 
ing boundary conditions: (i) At the surface r — R, we have 
the standard requirement AP — SP ~ pg^r = 0; (ii) At a 
small r close the the center, regularity requires SP/ p — Yqt'^ , 
= {f3/ujr)Yo r'^~^ {Yo is a constant), where /3 is deter- 
mined from f3{l + f3) = X, or (3 = (Vl + 4A- l) /2. Once 
LUr is found, the actual rotation rate is recovered from 

Sis = qUJrl'i'- 

Figure 2 shows the frequencies (in the inertial frame) of 
the j = m = 2 gl modes for model AU and UU as functions 
of the spin frequencies (These modes have Z = 2 in the f2s = 
limit). Note that as Q.s increases, ujr increases while LOi 
decreases. 

An approximate solution the the radial equation can 
be obtained by the WKB analysis. Substituting r^^r oc e''"^ 
and SPj p oc e^^"^ into the radial equations, we obtain the 
dispersion relation for g-modes: 

.2 A 



(11) 
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Figure 2. The frequencies of the j = m = 2 gi modes in the 
inertial frame for model AU and UU as functions of the spin 
frequency. The onset of instability occurs at i/; mode = 0. The 
solid lines are numerical results obtained by solving the radial 
equations; The dashed lines correspond to the WKB limits. 



Figure 3. Density perturbation Sp (in arbitrary unit) as a func- 
tion of radius for the j = m = 2 gi mode of model UU. The 
solid line is for Us = 0, and the dashed line is for Ps ~ 140 Hz 
(corresponding to q = 1.3). 



This indicates that ujr oc \/A. In fact, with the condition 
J^^ kdr = {n + C)tt (n is the number of radial node, C is 
a constant of order unity, ri and r2 are the inner and outer 
turning points where tj,- = N), we can write 

234Hz /Ay/y^^rflnA 

n + C \&) \ 300 Hz ) ^ ' 

Thus we have, for the j = 1 modes, a;, /'^''''' = -\/A/6 and 
fia/tj'''-' = q^J\/2A, where o;^"' is the corresponding mode 
frequency at zero rotation. The WKB results are also plot- 
ted in Fig. 2, showing remarkable agreement with the actual 
numerical result. Clearly, tj^ = occurs at 5 = 1, which 
corresponds to A = 11.13 from the solution of the m = 2 
Laplace equation (eq. |^). Thus, the onset of instability oc- 
curs when the ratio of spin frequency Q.s and the nonrotating 
mode frequency lo^^^ satisfies 

/ii.i3y/2 

The actual numerical calculations give slightly smaller criti- 
cal spin frequency (see Fig. 2). For Cls = Q.s/i^'"^^ < 1, in the 
WKB limit, the j = m = 2 mode frequency can be fitted to 
the following analytic expression to within 1%: ^ 

^ = + 1^^- 0.064 (14) 

Figure 3 gives two examples of the radial eigenfunctions. 
We see that the mode is concentrated in the outer region of 
the stellar core (0.6 ^ r/R <^ 0.9), and rotation tends to 

t A series solution (in q) of eq. ^ yields A = 6 + 2g + 10g^/7 + ■ ■ ■ 
for j = m = 2 (Ril dsten &; Ushomirsky 1996, private communi- 
cation). In eq. (|l4), the linear and quadratic terms are derived 
from this analytic expansion of A for small q and the last term is 
based on our numerical fitting. 



"squeeze" the mode into that region. This is because rota- 
tion increases ujr and narrows the region in which the wave 
can propagate [lOt < TV). The dispacements and ^± can 
be quite significant near the surface, but the associated en- 
ergy density (proportional to p|CP) is significant only in the 
propagation zone. 

We can similarly calculate the frequencies of higher- 
order g-modes. The results are even closer to the WKB 
limit (as expected). These modes tend to be unstable at 
lower spin frequencies. Similarly, the higher-m modes are 
also destablized by gravitational waves. We shall not be 
concerned with these higher-order modes since they have 
smaller growth rates due to gravitational radiation and 
higher damping rates due to viscosities as compared to the 
m = 2, gi mode (see Section 4). 

Finally, we discuss the validity of the two approxima- 
tions adopted in our calculations. As far as the mode fre- 
quency is concerned, the Cowling approximation is an ex- 
cellent approximation for g-modes: For nonrotating stars, 
our mode frequencies agree with those of Lai (1994) (where 
the full equations are solved) to within 2%. The Cowling ap- 
proximation does introduce an error (up to a factor of 80%) 
to the quadrupole moment of the mode (see §4). However, 
the uncertainty in the nuclear equation of state (particularly 
the nuclear symmetry energy which determines the g-mode 
properties) is presumably larger. As noted before, the con- 
ditions LUr <^ N and Qs ^ N, needed for the validity of 
the traditional approximation, are satisfied in the propaga- 
tion zone of the mode, although they are clearly violated 
near the stellar center and surface. Our numerical results 
indicate that for the gi mode, the inequality l^^] > is 
achieved only for r/R ^ 0.6 — 0.7 (depending on the rota- 
tion rate), and |^^| 3> \^r\ is achieved only near the stel- 
lar surface [indeed C±/Cr = GM/[R^ul) > 1 at r = i?]. 
Thus the traditional approximation is not rigorously justi- 
fied for the gi mode, although it becomes increasingly valid 
for higher-order modes. Clearly, without an exact, full cal- 



© 1999 RAS, MNRAS 000,|l|-| 



Secular Instability of g-Modes in Rotating Neutron Stars 5 



culation, which is currently beyond reach, it is difficult to 
know definitively how good the traditional approximation is 
when applying to the low-order modes. One way to assess 
the validity of our results is to consider the o;^"' ^ limit 
using the linear theory, which gives, for the I — m — 2 mode: 



J°'> - 2(1 - c„i)ns 



where 



C„2 



\dr 



(15) 



(16) 



(e.g., Unno et al. 1989 



^ l^rl, we find 



89) . If we assume l^i 
Cni = 1/6 and eq. (Ilq) reduces to the first two terms of 
eq. (^^. For model UU, our numerical calculation yields 
C„i = 0.104, 0.132, 0.148, 0.155 for gi,g2,g3,g4 respectively, 
corresponding to 2(1 - C„i) = 1-79, 1.74, 1.70, 1.69. These 
values are close to the corresponding value (5/3 = 1.67) 
obtained from the traditional approximation. Since our cal- 
culation indicates that the instability occurs at Q,s < tij'"-' 
(see eq. |l^), we are confident that our qualitative conclu- 
sion regarding the secular instability of g-modes is firm; the 
critical ratio fis/i^'"' for instability obtained in this paper 
(eq. |jl3[) is likely to be accurate to within ~ 20%. 



with 



£ = ^ I drr 



,SP{r) 



5p{r). 



(22) 



The dissipation due to gravitational waves is calculated 
from the multipole radiation formula derived by Lindblom 
et al. (1998) using the formalism of Thorne (1980): 



\ dt J , 



= -ujr}_^Niuji {\5DimV + \5Ji^\'). (23) 



l>2 

where A'^; is the coupling constant 
47rG {l + l){l + 2) 



Ni 



c2i+i 1(1 - 1)[(2/ + 1)!!]2 



(24) 



and the mass and current multipole moments of the per- 
turbed fluid are given by: 



5Dim = ld\5prYi*^, 



SJlm — - 



l + l 



d'^xr{piijJri + Spv) ■ Yi'^*, 



(25) 



(26) 



and is the magentic type vector spherical harmonic de- 
fined by Thorne (1980): 



4 DRIVING AND DAMPING OF THE MODES 

The condition that uji changes sign is only a necessary con- 
dition for the CFS instability. For the instability to be of 
interest, it must grow sufficiently fast and that the growth 
rate must overcome the possible damping rates due to vis- 
cosities. In this section we determine the mode growth rate 
and damping rate. 

The timescale r associated with any dissipative process 
is given by 

l = -±^ (17) 
T 2E dt' ^ ' 

where E is the canonical mode energy in the rotating frame 
(Friedman & Schutz 1978a): 



E 



+(v • c) r ■ (vp - c^vp) . 

We normalize the mode amplitude according to 
j d'xpt • C = ^ dr r^p {C + Aei) = 1, 



(18) 



(19) 



where A is evaluated using 

r-l 



A 



2tv 



dp 



1 (l^q^p^r 
m^(l+_aV) 2 



(l + qV)(l-M ) 



dHj„ 
dp 



AmqpHj 



\ dp 



(20) 



and the Hough function is normlized via 27r jdpHj^ = 1. 
The g-mode energy is dominated by the first term in eq. (|l8|), 
which gives 



E:^'^ul{l+£), 



(21) 



llm 



(27) 



For the j — m = 2 mode, the dominant gravitational 
radiation comes from the mass quadrupole SD22- We find 



25 



5 - -1 / SD22\ , ,3 „- 



(28) 



where cj = uj / [GM / B^Y^'^ , SD22 is dimensionless (in units 
where G = M = R = 1), Mia = M/(1.4Mq) and Ri = 
7?/(10km). Clearly, Tgr changes from positive (damping) for 
oji > to negative (driving) for uji < (where the mode is 
unstable for inviscid stars). 

Figure 4 shows the damping/growth rate Tg^.^ as a func- 
tion of the stellar spin frequency i^s for the gi mode of the 
two neutron star models. The change of with varying 
Us mainly comes from the changes in u>i and cOr, with the 
change in SD22 relatively minor. We see that when Vs ex- 
ceeds the critical spin rate for instability, the mode growth 
rate can be significantly larger than the damping rate for the 
nonrotating models (This is particularly evident for model 
AU). 

The higher order mass multipole radiations (from SD12, 
with / > 2) are less important. One can show that the ra- 
diation from 5J22 identically vanishes. The leading order 
current multipole radiation comes from 5J32, but it is much 
smaller than the mass quadrupole radiation. 

The growth/damping timescale of g-modes due to grav- 
itational radiation is typically much longer than that of f- 
modes (for which Tgr is of order seconds). This comes about 
for two reasons: (i) G-modes have lower frequencies; (ii) The 
multipole moment SDim for g-modes is much smaller than 
that for the f-modes (for which 5 Dim is of order unity). ^ 



3 The small SDim for g-modes can be understood from the sum 
rule (Reisenegger 1994) that must be satisfied for the multipole 
moments of all modes. 
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Figure 4. The gravitational radiation damping rate r^.^ as a 
function of spin frequency for the two g-modes depicted in Fig. 2. 
Positive Tg^^ (for smaU Us) imphes damping, while negative 
implies driving of the mode by gravitational radiation. At Us = 0, 
we have Tj 
for model AU. 



: 0.04 yr-i for model UU and t~'^ ~ 6 X lO'^ yi'^ 



The dissipation of mode energy due to viscosities can 
be calculated from (Ipser & Lindblom 1991) 



(-). 

\ dt / vise 



(29) 



where Saab and 5a are the shear and expansion of the per- 
turbation, respectively: 

2 



Saab = (^VaCb + Vb^a - -Sab^ cC^ 

5a — itiJrVc^'^. 



(30) 
(31) 



The shear viscosity resulting from neutron-neutron scatter- 
ing in the normal liquid core is given by (Flowers & Itoh 
1979; Cutler & Lindblom 1987): 



?7 = 1.1 X lO^Tg'^ 



9/4 

) g/(cm s), 

piiuc 



(32) 



where Tg is the temperature in units of 10 K, and pnuc = 
2.8 X 10^* g cm~^. The bulk viscosity results from the 
tendency to achieve chemical equilibrium via the modified 
URCA process during the perturbation (Sawyer 1989): 

2 



C = 4.8 X lO^^Tg^ 



P 

Pnuc 



UJr 

TP 



ll (cm s) 



(33) 



We have not been able to derive an explicit expression for 
[dE / dt)viBc from the g-mode wavefunctions to facilitate con- 
venient numerical computation of the viscous damping rate. 
We shall therefore be contented with order of magnitude es- 
timates. The mode damping rate due shear viscosity is ap- 
proximately given by the rate at which momentum diffuses 
across a mode wavelength: 

5/4 

yr"\ (34) 



'^shc 



Pnuc 



g-mode. Similarly, the mode damping rate due to bulk vis- 
cosity is estimated as 



Tbulk 



0. 



p 

Pnuc 



150 Hz 



(35) 



(where i/,- is the mode frequency in the rotating frame). The 
length scale of the mode is somewhat less than 10 km, as 
comparison with the numerical results for nonrotating g- 
modes (Lai 1994) indicates that eq. ( |34[ ) (with Li = 1) 
underestimates 1/rshear by a factor of ten, while eq. ( |35[ ) 
underestimates l/rbuik by a factor of three. When the star 
rotates faster than the mode frequency at zero rotation, the 
transverse wavelength of the mode decreases as ~ 2R/q (for 
g S> 1), and the radial wavelength also tends to be smaller 
(see Fig. 3). For the rotation rate of interest in this paper 
{Qs is at most a few times w'-^^), the mode wave function is 
only modestly changed by rotation, and we expect that rota- 
tion modifies the viscous damping rates by at most a factor 
of a few. Also note that the mode frequency (in the rotating 
frame) Vr increases with increasing spin frequency (e.g., for 
model UU, we have Ur = 296 Hz when i/s = 300 Hz), and 
this tends to reduce our estimate for t^^j^.. In any case, the 
uncertainty in our estimate of viscous damping rate is prob- 
ably comparable to the uncertainly in our understanding of 
the microscopic viscosity of neutron star matter. 

The total damping rate of the g-mode is given by 



1 _ 1 ^ 



1 1 

-I- 



''"shea 



Tbulk 



(36) 



< 0. Despite the uncer- 
ear and Tbulk, by comparing 



where L — lOLi km is the characteristic wavelength of the 



The mode is unstable when r 
tainty in our estimates for Tsh, 

Fig. 4 with eqs. ( ^^ and (^) we can readily conclude that 
for T ^ 10^ K the g-mode instability is suppressed by the 
bulk viscosity, while for T ^ 10^ K it is likely to be sup- 
pressed by the shear viscosity. For a narrow range of neutron 
star temperatures, around a few times 10^ K to 10® K (de- 
pending on the detail of viscous dissipation), the gi mode 
may be unstable due to gravitational radiation. To achieve 
this instability, the rotation rate of the star must be at least 
a factor of two to three greater than the mode frequency at 
zero rotation (see Fig. 4). 

So far we have only focused on the j = m = 2 gi mode. 
As the order of the g-modes increases, both the mode fre- 
quency and 5 Dim decreases, which tends to reduce iTg^-^]. 
Also, the viscous mode damping rate increases with increas- 
ing mode order. Therefore we expect that the CFS instabil- 
ities of higher-order g-modes are completely suppressed by 
viscous dissipations. 



5 VISCOSITY DRIVEN INSTABILITY? 

For an incompressible rotating Maclaurin spheroid, the I = 
—m = 2 f-mode becomes unstable due to viscosity at the 
same point as the CFS instability of the I — m — 2 mode 
(Chandrasekhar 1987). The I = —m > 2 viscous instabil- 
ity occurs at higher rotation rates. Recent studies indicate 
the instability occurs only for neutron stars with very stiff 
equation of state (Bonazzola et al. 1997). 

The viscous instability sets in when the mode frequency 
LOr in the rotating frame goes through zero. Figure 5 shows 
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Figure 5. The frequencies (in the rotating frame) of the j = 
—m = 2 gi modes for the two neutron star models (UU and AU) 
as functions of the spin frequency. For comparison, the dashed 
lines show the results for the j = —m = 3 gi modes. 



of changing the orbital decay rate (Reisenegger & Goldre- 
ich 1994; Lai 1994). Rotation may significantly change the 
strength of the resonance (Lai 1997), giving rise to measur- 
able effect in the gravitational wave form (Ho & Lai 1999). 
In another situation, when the neutron star is surrounded 
by a disk, tidal coupling between the stellar g-modes and 
disk density waves can drive the the unstable modes, which 
may affect the stellar rotation. By contrast, r-modes possess 
no mass quadrupole moment, and the coupling with disk is 
mediated only through higher order multipole moments and 
thus is expected to be weak. We wish to study some of these 
issues in a future paper. 
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the frequencies of the j = — m = 2 modes for the two neu- 
tron star models (UU and AU). For comparison, the fre- 
quencies of the j — —m = 3 modes are also shown. All 
these modes follow the approximate WKB scaling relation 
(see Sect. 3). We see that while LOr is reduced by 
rotation, it never goes to zero, up to the maximum possible 
rotation rate.|^ We therefore conclude that viscosity-driven 
instability is absent for g-modes. 



6 DISCUSSION 

We have shown in this paper that g-modes in an inviscid 
neutron star are susceptible to the CFS instability when the 
rotation rate is comparable to or greater than the mode fre- 
quencies. However, viscous dissipation tends to suppress the 
instability except possibly in a narrow range of temperatures 
around 10^ K. 

Comparing to f-modes, the g-mode instability has the 
advantage that it sets in at a lower rotation rate. Comparing 
to r-modes, the g-mode instability is much weaker. Calcula- 
tions (Lindblom et al. 1998; Andersson et al. 1998a) indicate 
that at spin frequency of 300 Hz, the growth rate of r-mode 
is approximately 10^ yr~^, more than three orders of magni- 
tude larger than 1 /rgr for the g-modes studied in this paper. 
Thus for isolated neutron stars, the r-mode instability plays 
an much more important role than the g-mode in deter- 
mining the spin evolution of the stars (Owen et al. 1998; 
Andersson et al. 1998b). 

However, one should not dismiss possible importance 
of neutron star g-modes and their stabilities in other situa- 
tions. In a merging neutron star binary, g-modes can be reso- 
nantly excited by the companion star, with the consequence 

^ Note that as the rotation rate approaches the maximum value, 
our calculation of the mode frequency breaks down since the ro- 
tational distortion becomes significant. 
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